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Abstract 

Certain  important  transcendental  equations  occur  in  the  case  of  elastic  angu¬ 
lar  regions  while  analyzing  these  regions  for  flexure,  vibration,  and  buckling. 
Previously  these  transcendental  equations  have  been  solved  for  roots  and  the 
data  has  been  tabulated  for  different  boundary  conditions,  as  the  angle  of  the 
region  is  varied.  The  purpose  of  this  paper  is  to  demonstrate  that  once  we 
solve  for  the  roots  at  a  speciflc  angular  region,  the  roots  for  angular  regions 
with  angles  ranging  from  0  to  27r  can  be  obtained  via  forward  integration. 

Key  Words:  Transcendental  equation,  elastic  angular  region,  flexure,  vibra¬ 
tion,  buckling  ,  polygonal  plate 
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1.  Introduction 

Certain  transcendental  equations  occur  in  the  case  of  solution  of  polygonal 
plate  problems  dealing  with  flexure,  vibration,  and  buckling  [1-3].  The  roots 
of  these  equations  need  to  be  evaluated  to  form  eigenfunctions  which  are 
essential  for  the  solution  of  the  mode  shape  of  the  plate  during  flexure  or 
vibration.  Also  the  eigenfunctions  can  be  used  to  find  the  frequencies  of 
vibration  of  sectorial,  triangular,  quadrilateral,  and  other  polygonal  plates. 

Several  different  methods  have  been  used  in  the  past  to  And  the  frequen¬ 
cies  of  vibration  of  plates.  An  overview  of  the  some  of  the  studies  on  the  free 
vibration  of  plates  is  given  in  [4].  Ref.  5  derives  an  equation  for  finding  the 
eigenfrequencies  of  polygonal  plates  with  free  simply  supported  mixed  edges 
and  Ref.  6  analyzes  the  free  vibration  of  right  triangular  plates  using  a  su¬ 
perposition  method.  In  [7],  the  Ritz  method  is  used  to  And  the  fundamental 
frequencies  of  five-sided  plates  which  are  obtained  by  cutting  out  an  isosceles 
triangle  from  one  corner.  The  dynamical  analogy  with  membranes  is  used 
in  [8]  to  study  the  free  vibration  of  regular  polygonal  plates  with  simply 
supported  edges.  In  [9]  finite  elements  are  used  to  analyze  annular  sectorial 
plates  having  their  inner  circular  edges  clamped.  All  of  these  plates  under 
arbitrary  boundary  conditions  can  be  analyzed  with  the  method  of  the  eigen¬ 
functions  given  in  [2].  Although  we  emphasize  the  classical  problems  related 
to  homogeneous  isotropic  thin  plates  under  small  deflections,  there  may  be 
suitable  extensions  of  the  methodology  to  the  anisotropic,  nonhomogeneous 
and  composite  material  cases  listed  in  [10]. 

Ref.  11  tabulates  the  first  ten  complex  roots  of  these  transcendental 
equations  as  the  angle  of  the  angular  region  is  varied  from  15°  to  180°  in  steps 
of  5°.  Also  if  real  roots  are  present,  these  are  also  tabulated  in  ascending 
order.  In  this  paper  we  show  that  once  the  roots  are  known  at  a  specific 
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angle,  the  roots  at  all  the  other  angles  ranging  from  0  to  27r  can  be  found 
by  solving  a  differential  equation  which  relates  the  roots  to  the  angle  of  the 
region. 


2.  Transcendental  Equations 

In  this  section,  we  indicate  how  the  transcendental  equation  for  an  angular 
region  arises  while  solving  for  basis  functions  for  the  clamped  case.  Once  the 
basis  functions  are  formed,  the  solution  for  flexure  or  vibration  of  the  angular 
region  can  be  expressed  as  a  linear  combination  of  the  basis  functions.  Since 
any  polygonal  plate  can  be  subdivided  into  angular  regions,  the  solution  on 
the  polygonal  plate  can  be  obtained  by  requiring  the  solutions  on  the  angular 
regions  to  satisfy  continuity  conditions  along  suitable  diagonals. 

Let  the  vertex  of  the  angular  region  be  at  (0, 0)  with  one  of  its  edges 
coinciding  with  the  x-axis.  Let  a  be  the  angle  of  the  angular  region.  The 
equation  of  motion  of  the  angular  plate  is  given  by  [12] 

p  d'^ip 


D  dt^ 


=  0 


where  D  is  the  flexural  rigidity,  p  is  the  uniform  mass  density  per  unit  area, 
V"  is  the  transverse  displacement,  and  =  V^V^,  being  the  Laplacian. 
Requiring  simple  harmonic  vibrations  V’  =  w{r,6)e~'^*,  we  get 


= 


Qj>2  Qj.  j.2  Q02 


W  =  pw 


where  p.  —  The  edges  of  the  angular  region  can  be  assigned  any  of  the 
classical  homogeneous  boundary  conditions.  As  an  example,  the  clamped 
boundary  conditions  are 

ti;(r,0)  =  0, 


w{r,  a)  =  0, 


dw 

'dd 


(r,  cc)  =  0. 
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For  the  sake  of  clarity,  we  now  indicate  how  basis  functions  and  transcen¬ 
dental  equations  can  be  formed  for  the  angular  region.  These  basis  functions 
are  used  in  [1]  for  finding  the  deflection  of  quadrilateral  plates  and  in  [2]  for 
finding  the  frequency  of  vibration  and  mode  shape  of  quadrilateral  plates. 

The  method  of  solution  for  =  fj,w  can  be  indicated  as  follows. 

Let  wo{r,0)  be  such  that  =  0,  and  wi{r,9)  be  such  that  =  /jlwq. 

Continuing  in  this  fashion,  we  define  Wi+i  by  =  ixwi.  Suppose  wi,  i  = 

1,2,...  satisfy  the  given  boundary  conditions.  Then  clearly  w  = 
satisfies  =  nw  and  the  boundary  conditions.  We  now  indicate  how  this 
procedure  can  be  applied  to  problems  associated  with  the  angular  region. 

We  have  =  0.  The  functions  wq  are  given  in  [1,11]  for  different 

boundary  conditions.  Note  that  on  any  compact  angular  region,  we  can  nor¬ 
malize  the  radius  vector  such  that  |r|  <  1  on  the  compact  region.  Choosing 

WQ^r,  O')  =  cos(A-t-l)^"t"5o  sin(A-|-l)0-f-C7o  cos(A  —  \)0-{-Dq  sin(A  1)^), 

the  clamped  boundary  conditions  imply  that  the  4x4  determinant 

/I  0  1  0 

cos(A4'l)a  sin(AH-l)Q!  cos(A— l)a  sin(A  — l)a 

det 

0  A-(-l  0  A-1 

\  — (A+1)  sm(A+l)a  (A+1)  cos(A+l)a  — (A— 1)  sm(A—l)a  (A  —  l)  cos(A— l)a 

The  above  equation  is  equivalent  to  the  transcendental  equation 

sin^  Aa  -  A^  sin^  a  =  0. 

Solving  the  transcendental  equation  numerically  for  A,  we  get  a  complex 
sequence  {Aj}^i  which  we  arrange  in  the  order  of  increasing  positive  real 
part. 

For  each  A,  we  can  determine  Ao,Bo,Co  and  Do  of  'woi'f',0)  upto  an 
arbitrary  multiplicative  constant  (e.g.,  let  Aq  =  1)  based  on  the  boundary 
conditions. 
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Let  A  =  Aj  for  some  i.  We  need  to  get  w\  such  that 

=  IJ.WQ  =  cos(A  +  1)0  +  Bq  sin(A  +  1)0 

+C'ocos(A  -  1)0  +  I>osin(A  -  1)0). 


Hence, 


Aq  cos(A  +  1)0  +  Bq  sin(A  +  1)0 
[(A  +  5)2-(A+1)21[(A  +  3)2-(A  +  1)21 
Cocos(A  —  1)0  +  £)osin(A  —  1)0 
[(A  +  5)2-(A  +  1)2|[(A  +  3)2-(A  +  1)2]. 

cos(A  +  5)0  +  Bi  sin(A  +  5)0 
+Ci  cos(A  +  3)0  +  Di  sin(A  +  3)0], 


where  the  homogeneous  part  cos(A+5)0+5i  sin(A  +  5)0+Ci  cos(A  + 

3)0  +  I?i  sin(A  +  3)0]  enables  satisfaction  of  the  boundary  conditions  on  0  =  0 
and  0  =  a. 

Thus  select  Ai,Bi,Ci,  and  £»i  to  satisfy  the  four  boundary  conditions 

u;i(r,0)  =  u;i(r,Q:)  =  ^^(’’.0)  = 

Since  A  satisfies  the  transcendental  equation  sin^  Acc  —  A^  sin^  a  =  0,  to  get 
a  nonsingular  system  for  the  solution  of  Ai,Bi,Ci,  and  Di,  we  should  have 
sin^(A  +  4)a  -  (A  +  4)^  sin^  a  ^0.  This  can  be  easily  verified. 

Continuing  in  this  manner,  we  get  for  i  >  1, 

Wi  =  wi  +  cos(A  +  1  +  4i)0  +  Bi  sin(A  +  1  +  4i)0 

+Ci  cos(A  -  1  +  4i)0  +  Di  sin(A  -  1  +  4i)0), 


where  V'^Wi  =  fiwi-i,  and  Ai,Bi,Ci,  and  Di  are  selected  such  that  Wi  sat¬ 
isfies  the  boundary  conditions  along  the  edges. 
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3.  Differential  Equations 

Let  a  be  the  angle  of  the  angular  region  and  u  be  the  Poisson’s  ratio.  The 
transcendental  equations  for  various  boundary  conditions  are  given  by  the 
following  table  [1,11]- 


TABLE  1:  Transcendental  Equations 

Boundary  Condition 

Case 

9  =  0 

9  =  a 

Transcendental  Equation 

1 

Clamped 

Clamped 

sin^  Aa  =  A^  sin^  a 

2 

Free 

Free 

(3  +  v)^  sin^  Aa  =  (1  —  v)^X^  sin^  a 

3 

Clamped 

Free 

(3  -f  i/)(l  —  i/)sin^  Aa  =  4  —  (1  —  i/)^A^sin^a 

4 

Clamped 

SS 

sin  2Aa  =  A  sin  2a 

5 

SS 

Free 

(3  -t-  v)  sin2Aa  =  —(1  —  u)X  sin  2a 

6 

SS 

SS 

sm  Aa  =  sm  a 

In  the  table  above,  SS  means  simply  supported.  In  cases  1,  2  and  6, 
the  transcendental  equations  may  be  further  broken  down  into  those  for 
symmetrical  and  anti-symmetrical  modes. 

In  cases  1,2,  and  6,  the  transcendental  equation  can  be  written  as 

sin^  \oL  =  {bX  +  c)^  sin^  a 

for  suitable  constants  a,  b  and  c.  The  transcendental  equation  for  the  sym¬ 
metrical  mode  is 

a  sin  Xa  +  {bX  +  c)  sin  a  =  0. 

Taking  differentials  on  both  sides,  we  get 

[A  da  +  a  dX]a  cos  Xa  +  bdX  sin  a  -f-  {bX  +  c)  cos  ada  —  0. 
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This  results  in 

dX  —  [aAcosAo!  +  (6A +  c)cosa] 
da  aa  cos  Xa  +  b  sin  a 

In  the  anti-symmetrical  case,  we  get 

dX  _  -  [aA  cos  Aa  -  (6A  +  c)  cos  a] 
da  aa  cos  Xa  —  b  sin  a 

In  cases  4  and  5,  the  transcendental  equation  can  be  put  in  the  form 

a  sin  2Aq;  =  bX  sin  2a 

for  suitable  a  and  b.  The  differential  equation  can  be  derived  as 

dX  2(6Acos2a  —  aAcos2AQ:) 
da  2aa  cos  2Aa'  —  b  sin  2a 

For  case  3,  the  transcendental  equation  is 

a  sin^  Aa  =  4  -  bX^  sin^  a, 

where 

a  =  (3  +  v){\  -  u), 

6  =  (1-W". 

The  relevant  differential  equation  is 

dX  —  [6A^sin2a  +  aAsin2Aa] 
da  aa  sin  2Aa  +  2bX  sin^  a 

5 

The  parameter  ^  —  Xa  tends  to  vary  more  slowly  than  A  as  a  is  varied 
from  0  to  27r.  Thus,  from  a  computational  point  of  view,  it  might  be  beneficial 
to  formulate  the  differential  equations  in  terms  of  /i.  For  the  symmetrical 
mode  of  cases  1,2,  and  6,  the  differential  equation  is 

dji  —  (6/i -|- ca)  cos  a -h  ^  sin  a 
da  aa  cos  /x  -f-  6  sin  a 
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For  the  anti-symmetrical  mode  of  cases  1,2,  and  6,  the  diflFerential  equation 
is 

dfji  (bfi  +  ca)  cosa  —  ^  sin  a 
da  aa  cos  fi  —  b  sin  a 

The  differential  equation  for  cases  4  and  5  is 

dfj,  b{2fj,  cos  2a  -  ^  sin  2a) 
da  2aa  cos  2/x  —  6  sin  2a 

For  case  3,  the  differential  equation  is 

d/x  asin2a -h  2sin^a] 

da  aa  sin  2 fj,  +  2b  sin^  a 

It  can  be  observed  from  the  numerators  of  the  right  sides  of  the  differen¬ 
tial  equations  for  /x  that  as  a  tends  to  zero,  the  numerators  in  all  cases  tend 
to  zero  also.  Thus,  for  moderate  a,  the  parameter  /x  is  relatively  constant. 

4.  Advantages  of  Forward  Integration 

One  of  the  advantages  of  the  proposed  method  is  that  if  the  roots  are  found  at 
an  arbitrary  angle,  the  roots  at  any  other  angle  in  the  range  of  the  initial  and 
final  angles  can  also  be  tabulated  with  significantly  reduced  computing  effort. 
Only  forward  integration  is  needed  to  compute  the  corresponding  roots  at 
the  intermediate  angle.  However,  as  can  be  observed  from  the  examples  in 
Section  5,  the  forward  integration  may  not  yield  all  the  roots  and  in  fact  may 
occasionally  result  in  spurious  roots.  We  observed  that  almost  all  of  these 
missing  roots  are  real  and  it  is  not  difficult  to  compute  these  using  any  of 
the  standard  methods.  Some  extra  effort  is  needed  to  compute  the  missing 
complex  roots  and  these  can  be  computed  using  a  search  procedure.  In  spite 
of  these  obstacles,  there  is  a  significant  reduction  in  the  computing  effort  by 
the  use  of  the  forward  integration  technique  since  the  usual  search  procedure 
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for  all  the  complex  roots  is  more  time  consuming  and  involves  an  intelligent 
choice  for  the  initial  approximation  for  the  roots. 

For  moderate  a,  since  the  parameter  yu  =  Aa  is  relatively  independent  of 
a,  the  time  taken  for  the  initial  approximation  of  the  roots  can  be  reduced. 
However,  when  we  consider  the  whole  range  of  [0, 27r],  the  value  of  fi  varies 
significantly.  We  present  several  examples  in  Section  5  illustrating  the  vari¬ 
ation  of  /X  with  a.  Our  method  is  especially  useful  in  cases  where  there  is 
significant  variation  in  n  with  a..  It  is  easy  to  plot  the  various  values  of  /x 
as  a  varies. 

Almost  all  the  complex  roots  and  most  of  the  real  roots  are  obtained 
by  this  procedure.  However,  some  of  the  roots  need  to  be  found.  It  is  a  lot 
simpler  to  find  the  real  roots.  One  of  the  observations  from  Ref.  11  is  that 
at  least  for  0  <  a  <  tt,  the  real  roots  are  smaller  than  the  least  magnitude 
of  the  real  part  of  the  complex  roots  for  almost  all  types  of  boundary  condi¬ 
tions.  The  only  exception  to  this  is  the  clamped-free  case.  We  will  present 
examples  involving  clamped,  clamped-free  and  simply  supported-free  bound¬ 
ary  conditions  in  Section  5  to  illustrate  the  usefulness  of  our  approach.  It 
happens  on  occasion  that  the  forward  integration  generates  spurious  roots. 
Thus  it  is  beneficial  to  check  the  validity  of  the  roots  obtained  by  the  forward 
integration  procedure. 

Another  significant  advantage  with  the  approach  of  the  differential  equa¬ 
tions  is  that  the  initial  value  problem  can  be  integrated  over  the  range  of 
a  from  e  to  27r,  where  e  approaches  0.  The  differential  equations  become 
singular  at  a  =  0.  The  roots  as  a  approaches  27r  are  useful  for  the  solution 
of  the  crack  problem. 
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5.  Examples 

We  solve  three  examples  in  this  section  to  demonstrate  the  usefulness  of 
the  approach  in  the  case  of  the  three  different  types  of  differential  equations 
considered  in  Section  3. 

Example  1:  Let  us  take  the  transcendental  equation  for  the  symmetrical 
roots  given  by 

sinAa  +  A  sin  a  =  0 

for  an  angular  region  with  clamped  boundary  conditions.  Note  that  p  =  Xa. 
For  a  =  15°,  the  first  sixteen  roots  of  the  transcendental  equation  in  the 
order  of  increasing  positive  real  part  and  the  corresponding  values  of  ^  are 
evaluated  and  are  given  in  Table  2. 

Let  p,  =  p/tt.  To  get  the  value  of  ^  at  a  =  60°,  the  differential  equation 

dp  —pcosa  +  ^sina 
da  a  C0S7r/i  +  since 

was  solved  with  the  initial  condition  /i(7r/12)  given  by  the  values  at  a  = 
15°  given  above.  We  used  the  fourth  order  Runge-Kutta  method  given  in 
the  MATLAB  package  using  a  tolerance  of  le  —  10.  The  values  of  /x/tt  at 
a  =  60°,  175°,  225°  and  355°  obtained  by  the  forward  integration  method 
are  listed  in  Table  3.  The  real  roots  that  were  not  generated  by  the  forward 


integration  are  listed  in  Table  4. 
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TABLE  2:  Values  of  A  and  n/ir  a.t  a  =  15°  for  the  Clamped  case 


A 

16.099036079369+  8.549629866614  i 
40.922739395410  +  11.808479416584  i 
65.217933871100  +  13.519976459038  i 
89.376965549419  +  14.695568014823* 
113.478110132424  +  15.592891970920* 
137.548757675389  +  16.318980707314* 
161.601199318519  +  16.928883919527* 
185.641832424667  +  17.454722092831  * 
209.674337397439  +  17.916891466163* 
233.700991180886  +  18.329159178917* 
257.723282631706  +  18.701258444435* 
281.742228750921  +  19.040325717266* 
305.758549248677  +  19.351749656939* 
329.772768581551  +  19.639699129280* 
353.785278459033  +  19.907465712142* 
377.796377663910  +  20.157693813903* 


_ /fA _ 

1.341586339947  +  0.712469155551  * 
3.410228282951  +  0.984039951382  * 
5.434827822592  +  1.126664704920* 
7.448080462452  +  1.224630667902* 
9.456509177702  +  1.299407664243* 
11.462396472949  +  1.359915058943* 
13.466766609877  +  1.410740326627* 
15.470152702056  +  1.454560174403* 
17.472861449787+  1.493074288847* 
19.475082598407  +  1.527429931576  * 
21.476940219309  +  1.558438203703* 
23.478519062577  +  1.586693809772  * 
25.479879104056  +  1.612645804745* 
27.481064048463  + 1.636641594107* 
29.482106538253  +  1.658955476012* 
31.483031471993  +  1.679807817825  * 


The  values  of  /*  in  Table  3  were  verified  by  solving  the  transcendental 
equation  sin  A  a  +  A  sin  a  =  0  for  a  =  60°  and  the  resulting  match  was  at 
least  to  the  tenth  decimal  place. 

Example  2:  For  the  clamped-free  case,  with  v  =  0.3,  the  transcendental 
equation  is 

2.31  sin^  Aa  =  4  —  0.09A^  sin^  a. 

The  values  of  A  and  p/tt  for  a  =  15°  are  listed  in  Table  5.  The  differential 
equation 

dp  0.497r(^)^  [— a  sin2Q:  +  2sin^  a] 

2.31asin27r)U  +  0.98(^)sin2a 

was  solved  with  the  initial  condition  p{'k/ 12)  and  the  resulting  values  of  p  for 
a  =  60°  and  a  =  160°  are  tabulated  in  Table  6.  Table  7  lists  the  additional 
roots  that  were  missed  by  the  forward  integration  procedure.  There  was  also 
a  spurious  value  of  /*  =  1.550263  +  0.000000*  generated  at  a  =  160°  that 
does  not  correspond  to  a  root  of  the  relevant  transcendental  equation. 


TABLE  3:  Values  oi  ii/it  for  various  angles:  Clamped-Claraped  (Sym.)  Case 

a  =  60°  I  a  =  175°  I  a  =  225°  I  a  =  355° 

1.353109  +  0.650683  i  1.944236  +  0.000000  i  2.449112  +  0.277056  i  2.028574  +  0.000000  i 

3.415242  +  0.925988  i  3.886716  +  0.000000  i  4.461036  +  0.504582  i  4.057384  +  0.000000  i 
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9.314538  +  0.000000  i  8.871763  +  0.000000  i 

10.840970  + 0.000000  i 
12.808479  + 0.000000  i 
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TABLE  5:  Values  of  A  and  /^/tt  at  a  =  15°  for  the  Clamped-Free  case 

A  p/'n 

10.651024802745+  0.000000000000  i  0.887585400229  +  0.000000000000  i 

7.610592426966+  1.501062004362  i  0.634216035580  +  0.125088500364  z 

22.784846703345+  6.342239495230  z  1.898737225279  +  0.528519957936  z 

35.073902145999  +  8.093194426143  z  2.922825178833  +  0.674432868845  z 

47.233021440439  +  9.253448430343  z  3.936085120037  +  0.771120702529  z 

59.338504245872  +  10.131592314675  z  4.944875353823  +  0.844299359556  z 

71.414977259846  +  10.840846669030  z  5.951248104987  +  0.903403889086  z 

83.473543259896  +  11.436757789937 z  6.956128604991  +  0.953063149161 z 
95.520120582280  +  11.951029480485  z  7.960010048523  +  0.995919123374  z 
107.558209383922  +  12.403569139462 z  8.963184115327  +  1.033630761622  z 
119.590034126021  +  12.807735494723 z  9.965836177168  +  1.067311291227 z 
131.617085967549  +  13. 172944819256  z  10.968090497296  +  1.097745401605  z 
143.640406736456  +  13.506091735206  z  11.970033894705  +  1.125507644601  z 
155.660748892999  +  13.812380105128  z  12.971729074417  +  1.151031675427  z 
167.678670940916  +  14.095836412439  z  13.973222578401  +  1.174653034370  z 
179.694597021781  +  14.359641355224  z  14.974549751815  +  1.196636779602  z 
191.708855583596  +  14.606352014275  z  15.975737965210  +  1.217196001190  z 


TABLE  6:  Values  of  /z/zr  f 

or  the  Clamped-Pree  Case 

P 

II 

o 

o 

a  =  160° 

1.079931  + 0.000000  z 

1.744045  (false  root) 

0.567564  + 0.190356  z 

0.501101  + 0.245940  z 

1.904356  + 0.457347  z 

2.505965  + 0.216165  z 

2.927638 +  0.611914Z 

3.509140  + 0.178136  z 

3.940002  + 0.711132  z 

4.513539  +  0.096257Z 

4.948141  + 0.785436  z 

5.363897  + 0.000000  z 

5.954037  + 0.845148  z 

6.251376  + 0.000000  z 

6.958558  + 0.895174  z 

7.121553  + 0.000000  z 

7.962160  + 0.938271  z 

7.978357  +  0. 159757  z 

8.965111  + 0.976149  z 

8.980855  + 0.257304  z 

9.967581  +  1.009949  z 

9.982578  + 0.320489  z 

10.969685  +  1.040473  z 

10.983871  + 0.369089  z 

11.971501 +  1.068304Z 

11.984896  + 0.409186  z 

12.973088  +  1.093882  z 

12.985740  + 0.443597  z 

13.974487  + 1.117546  z 

13.986454  + 0.473889  z 

14.975733  +  1.139565  z 

14.987071  + 0.501033  z 

15.976849  +  1. 160153  z 

15.987613  + 0.525682  z 
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Example  3:  For  the  simply  supported-free  case,  with  v  =  0.3,  the  transcen¬ 
dental  equation  is 

3.3sin2AQ:  =  — 0.7Asin2Q:. 

The  values  of  A  and  p/i:  for  a  =  15°  are  listed  in  Table  8. 


TABLE  8:  Values  of  A  and  /^/tt  at  a  =  15°  for  the  SS-Free  case 

_ A _ /x/tt _ 

7.895591872816  -I-  0.000000000000  i  0.657965989401  0.000000000000  i 

9.307012034278  +  0.000000000000 i  0.775584336190  +  0.000000000000 i 
20.718871978887  -I-  2.743152513018  i  1.726572664907  +  0.228596042752  i 
32.777489190922  +  3.676260674920 i  2.731457432577 -b  0.306355056243 i 
44.813676592684  -b  4.288658065852 i  3.734473049390  -b  0.357388172154 i 
56.838604582663  +  4.748467761566 i  4.736550381889  -b  0.395705646797 i 
68.856982193617  +  5.117556866893 i  5.738081849468  -b  0.426463072241  i 
80.871175036583  -b  5.426171354656 i  6.739264586382  -b  0.452180946221 i 
92.882514190978  -b  5.691485721814 i  7.740209515915  +  0.474290476818 i 
104.891811137801  -b  5.924226510378 i  8.740984261483  -b  0.493685542532 i 
116.899591164329  +  6.131555242723 i  9.741632597027 -b  0.510962936894 i 
128.906210549635  +  6.318496510883 i  10.742184212470  +  0.526541375907 i 
140.911920160849  -b  6.488713926984*  11.742660013404  -b  0.540726160582* 

152.916902153721  -b  6.644962269170  *  12.743075179477  -b  0.553746855764  * 

164.921292234989  -b  6.789365833965  *  13.743441019582  -b  0.565780486164 * 

176.925193736710  -b  6.923597589457*  14.743766144726  -b  0.576966465788* 

188.928686829963  -b  7.048998802756  *  15.744057235830  -b  0.587416566896  * 

200.931834739295  +  7.166661433097*  16.744319561608  -b  0.597221786091  * 

212.934688043977  +  7.277486401963  *  17.744557336998 -b  0.606457200164  * 

224.937287723129  +  7.382225757249  *  18.744773976927  -b  0.615185479771  * 
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TABLE  9:  Values  of  ju/7r  for  the  SS-Pree  Case 

a  =  60° 

a  =  160° 

0.548912  + 0.000000  i 

0.488073  +  0.000000* 

0.915801  + 0.000000  i 

1.025131  +  0.000000* 

1.780755  +  0.000000 i 

2.050935  +  0.000000* 

2.738015  + 0.155043  i 

3.078267  +  0.000000* 

3.739534  + 0.215478  i 

4.108517  +  0.000000* 

4.740664  + 0.257281  i 

5.144743  +  0.000000* 

5.741543  + 0.289754  i 

6.199987  +  0.000000* 

6.742251  + 0.316452  i 

7.246253  +  0.074584* 

7.742834  +  0.339177* 

8.246425  +  0.113503* 

8.743324  +  0.358984* 

9.246579  +  0.140984* 

9.743743  +  0.376552* 

10.246719  +  0.163023* 

10.744106  +  0.392343* 

11.246845  +  0.181689* 

11.744424  +  0.406688* 

12.246961  +  0.197998* 

12.744705  +  0.419833* 

13.247067  +  0.212541* 

13.744956  +  0.431965* 

14.247165  +  0.225700* 

14.745181  +  0.443229* 

15.247255  +  0.237739* 

15.745384  +  0.453746* 

16.247339  +  0.248847*. 

16.745569  +  0.463602* 

17.247417  +  0.259168* 

17.745738  +  0.472882* 

18.247490  +  0.268812* 

18.745892  +  0.481647* 

19.247558  +  0.277869* 

TABLE  10:  Additional  values  oi  p/n:  SS-Pree  case 

a  =  60° 

a  =  160° 

5.072741  +  0.000000* 

1.463953  +  0.000000* 
2.050935  +  0.000000* 
3.412327  +  0.000000* 
4.382660  +  0.000000* 
6.292132  +  0.000000* 

The  differential  equation 

dp  0.7 [sin  2a  —  2a  cos  2a] 
da  6.6a  cos  27r/i  +  0.7  sin  2a 

was  solved  with  the  initial  condition  /i(7r/12)  and  the  resulting  values  of  p  for 
a  =  60°  and  a  =  160°  are  tabulated  in  Table  9.  Table  10  lists  the  additional 
roots  that  were  missed  by  the  forward  integration  procedure. 
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6.  Conclusions 

The  roots  of  the  transcendental  equations  can  be  found  for  any  arbitrary 
angular  region  by  solving  the  associated  initial  value  problem.  This  proce¬ 
dure  is  less  time  consuming  since  the  alternate  method  of  solving  the  original 
transcendental  equation  requires  a  search  procedure  in  the  neighborhood  of 
the  roots. 
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